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ABSTRACT 

We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for 
low-mass compact binary coalescence (CBC) searches during LIGO's sixth science run and Virgo's 
second and third science runs. We present strain noise power spectral densities (PSDs) which are 
representative of the typical performance achieved by the detectors in these science runs. The data 
presented here and in the accompanying web-accessible data files p] are intended to be released to 
the public as a summary of detector performance for low-mass CBC searches during S6 and VSR2-3. 

I. INTRODUCTION 

The LIGO and Virgo gravitational wave observatories recently completed their final joint science run until advanced 
detectors come online. The observation was LIGO's sixth science run and Virgo's second and third science runs. LIGO 
operated two instruments, the four kilometer Hanford detector (HI) and the four kilometer Livingston detector (LI) 



from July 07, 2009 to October 20, 2010. The two kilometer Hanford interferometer, which was included in the 
S5/VSR1 CBC analyses, was decommissioned prior to the start of this run and did not collect any data. The three 
kilometer Virgo interferometer (VI) operated during roughly the same timespan, with however a major commissioning 
break between its second and third science runs, from Jan 11, 2010 to August 6, 2010. 

In this note, we summarize the sensitivity achieved by these detectors during the latest runs from the perspective of 
low-mass (inspiral-only) CBC searches. The strain noise power spectral density (PSD) is a complete characterization 
of the sensitivity of a detector but generally only meaningful over timescales of about an hour. Over longer time scales, 
the noise characteristics of the detectors typically vary significantly, for instance as the morning traffic picks up. More 
importantly, the detectors themselves change over time as they undergo weekly maintenance and occasional larger 
scale upgrades. The notion of PSD does not immediately have meaning when applied to the detector performance 
over the entire duration of the run. To solve this problem, we choose a single "representative" block of time for each 
detector and compute the PSD for each detector in this block. We take the resulting PSD as representative of the 
typical sensitivity to gravitational waves from CBCs achieved by the detectors in S6/VSR2-3. 

We choose the respresentative time for each detector by looking at its inspiral horizon distance distribution and 
selecting a time for which the detector operated with a horizon distance close the the mode of this distribution. 
The inspiral horizon distance is equal to the largest distance at which an equal-mass compact binary inspiral would 
accumulate an average SNR of 8 in the detector. The inspiral horizon distance contains less information but is derived 
from the PSD and includes information specific to CBC signals. In this sense, the horizon distance is a useful measure 
of the sensitivity of a detector to gravitational waves from CBCs at a given time. In this article, we gather the inspiral 
horizon distance data generated during S6/VSR2-3 low-mass CBC search [2] and use the results to identify PSDs that 
are representative of detector performance for CBC searches during these science runs. 

The plots and data presented here are intended to be released to the public as a summary of detector performance 
for CBC searches during S6 and VSR2-3 [Tj. These results use exactly the same science segments and analysis code 
[3] that was used in the low-mass CBC search in S6 and VSR2-3. In the next section, we review how the inspiral 
horizon distance is computed from the PSD. In section 3, we present the inspiral horizon distance data obtained 
from S6/VSR2-3 CBC analyses. We then use these data to identify PSDs which are representative of the detector 
sensitivity to inspirals for this science run. We close this note with a discussion of one potential pitfall in using these 
spectra. 

II. INSPIRAL HORIZON DISTANCE 

In LIGO and Virgo data analysis applications, we treat the strain noise in a detector as a stationary random process. 
If the noise in the detector were truly stationary, then the noise spectral density would completely characterize the 
sensitivity of the detector as a function of frequency. The power spectral density S n (f) for a stationary random 
process n(t) is defined implicitly by the relation 

ls n (f)S(f-f') = {n(f)n*(f)), (1) 

where h(f) is the Fourier transform of the random process. The spectral density is a measure of the mean square noise 
fluctuations at a given frequency. As mentioned above, the noise in the LIGO and Virgo detectors is not stationary. 
However, by measuring the spectral density over a short enough timescale, we are able to approximate the noise as 
stationary. The chosen timescale must also be long enough that we can form an accurate estimate of the spectral 
density. In the S6/VSR2-3 CBC searches, the spectral density was computed on 2048-second blocks of contiguous 
data |2]. We account for long timescale non-stationarities by using a different spectral density for every 2048 seconds. 
In assessing the overall performance of a detector for CBC searches, we use the inspiral horizon distance data from 
S6 and VSR2-3 to identify the "typical" sensitivity of the interferometers. The inspiral horizon distance of a detector 
is the distance at which an optimally oriented and optimally located equal-mass compact binary inspiral would give 
an average signal to noise ratio (SNR) of p = 8 in the interferometer. If h(f) represents the Fourier transform of the 
expected signal, then the average SNR this signal would attain in a detector with spectral density S n (f) is given by 




We find the inspiral horizon distance by setting (p) — 8 and solving for the distance D to the inspiral event which 
parametrizes the waveform h(f). Thus, the inspiral horizon distance combines the spectral density curve with the 
expected inspiral waveform to produce a single quantity that summarizes the sensitivity of the detector at a given 
time. 



Practical considerations require modifications to the limits of the integral. In the CBC search code, we compute 
the signal to noise ratio by 




rA..» \h(fW 

<P> = l/4 L ~sM) df - (3) 

The lower limit is determined by our ability to characterize the noise at low frequencies. In the S6 CBC search, we 
took f[ ow = 40Hz as the low frequency cut-off in computing the inspiral horizon distance. For Virgo in VSR2-3, the 
low frequency cut-off was fi ow = 50Hz. The upper limit of the integral is the innermost stable circular orbit (ISCO) 
frequency, 

c 3 
flsc ° = GV&rGM' (4) 

where M is the total mass of the binary system. For binary neutron star (BNS) systems, for which we take m = 1.4M©, 
f isco = 1570Hz. 

The inspiral waveform for CBCs is accurately given in the frequency domain by the stationary phase approximation. 
For an optimally oriented and optimally located equal mass binary, the signal that appears at the interferometer (in 
this approximation) is given by 

HI) = ± (^3) V2 (GMf/ e (nf)-y^f^\ (5) 

where M. = fj?/ 5 M 2 / 5 is the chirp mass of the binary, D is the distance to the binary and W is a real function of /, 
parametrized by the total mass M. Setting (p) = 8 and inserting this waveform into eqn. k3l we find that the inspiral 
horizon distance is given by 



*=K£f ^'WC'G^- (6) 

The inspiral horizon distance is defined for optimally located and oriented sources. To compare with previous results, 
note that the sensitive range of an interferometer gravitational-wave detector was considered by Finn and Chernoff 
[5], though here we following the conventions of Allen, et. al. |4j and Brown [6]. Furthermore, if we divide the inspiral 
horizon distance given here by 2.26 we obtain the SenseMon range 7] reported as a figure of merit in the LIGO and 
Virgo control rooms, where the factor of 2.26 comes from averaging over a uniform distribution of source sky locations 
and orientations. 

In practice, it is convenient to measure distances in Mpc and mass in Mq. It is useful therefore to specialize eqn. 
m\ to this unit system. Further, since we measure the strain h(t) at discrete time intervals At = l// s , the spectral 
density is only known with a frequency resolution of A/ = f s /N, where N is the number of data points used to 
measure S n (f). By putting / = k/(NAt) into eqn. pland grouping terms by units, we arrive at the expression 
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for the inspiral horizon distance in Mpc. Since it is convenient to work with the binary system's component masses, 
we have also replaced the chirp mass M. with the reduced mass \i and the total mass M. Written this way, the inspiral 
horizon distance in Mpc is easily computed from the binary component masses in Mq. 

III. SUMMARY OF INSPIRAL HORIZON DISTANCE DATA 

Here we present the horizon distance data collected from the final data products produced during the S6/VSR2-3 
lowmass CBC search [2]. We have collected the data, rather than computing the inspiral horizon distance directly, in 
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FIG. 1: (a) Inspiral horizon distance as a function of time during S6-VSR2/3. The average inspiral horizon distances for each 
week in S6 and VSR2-3. As an indication of the weekly variations, we have included error bars corresponding to the standard 
deviation of the inspiral horizon distance during each week, (b) Distribution of 1.4-1.4 solar mass inspiral horizon distance for 
the three gravitational wave detectors HI, LI, and VI for the joint LIGO- Virgo science run consisting of S6 and VSR2/3. The 
histogrammed data consists of the same 2048-second analyzed segments from the S6 and VSR2/3 CBC searches. 

order to ensure that we analyze the exact same science segments and use the exact same analysis code as used in the 
LIGO/Virgo CBC searches. 

In Fig. [ij,, we plot the average BNS inspiral horizon distance for each of the three detectors as a function of time. 
We use a window of one week and the points on the plot correspond to the average inspiral horizon distance for all 
science segments beginning in that week. The error bars attached to the points indicate the standard deviation in the 
inspiral horizon over the course of the given week. This figure highlights the variability in sensitivity throughout the 
run and the reason it is difficult to identify a single time for each detector with a typical or average sensitivity. In Fig. 
[TJd, we histogram the BNS inspiral horizon distance for the three detectors HI, LI, and VI. The bimodal behavior 
seen in the LIGO and Virgo detectors is largely due to a significant commissioning break in S6 and commissioning in 
Virgo between VSR2 and VSR3. These commissioning breaks will be described in detail in a later publication on the 
S6/VSR2-3 runs. 

In the actual S6/VSR2-3 CBC analysis, the inspiral horizon is computed for (n)-(n) solar mass binaries for n an 
integer. Previous documents [5] however have plotted the horizon distance for the canonical 1.4— 1.4 solar mass binary 
neutron star. In order to simplify comparison to previous results, we rescale the obtained distributions by (2.8/2) 5 / 6 
corresponding to the ratio of chirp masses of a 1.0 — 1.0 solar mass system and a 1.4 — 1.4 solar mass system. This 
scaling ignores the fact that fi SCO is different for the two mass pairs, but this is negligible since the signal template is 
buried in the noise at such high frequencies. 

In Fig. [2] we show the mean inspiral horizon distance for each interferometer as a function of the binary total 
mass, assuming equal mass binaries. This plot reflects the mean performance of the detector over various frequency 
bands. As the component mass becomes higher, the upper cutoff frequency fhigh — fisco becomes smaller and smaller. 
This means that the inspiral horizon distance focuses on a narrower band around the lower cutoff fi ow = 40Hz (or 
flow — 50Hz in the case of Virgo). The inspiral horizon distance takes into account only the inspiral stage of the 
CBC event, while for high-mass systems (M > 25M Q ) the merger and ringdown stages of the coalescence occur in 
the LIGO and Virgo sensitive band. For total masses greater than 25M Q , the inspiral-only range begins to fall over, 
which is not indicative of the sensitivity of the detector for these systems. For these binary systems, we use Effective 
One Body Numerical Relativity (EOBNR) waveform templates that include the merger and ringdown stages and our 
sensitivity is significantly improved relative to an inspiral-only analysis [5]. 



IV. REPRESENTATIVE POWER SPECTRAL DENSITY 



In Fig. [3j we give representative spectral density curves for each of the three detectors during S6 and VSR2-3. The 
chosen representative curve corresponds to a time when the detector operated near the mode of its inspiral horizon 
distance distribution shown in Fig. [TJd. The algorithm used to compute the spectral densities is described in detail 
in [5] . The parameters needed in order to reconstruct our results are given in Table [[] The first column in Table 
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FIG. 2: Mean inspiral horizon distance as a function of mass for the three gravitational wave detectors HI, Ll and VI during 
S6-VSR2/3. The error bars on the curves extend from one standard deviation below to one standard deviation above the mean. 



TABLE I: Parameters used in the computation of the spectral density. 



FINDCHIRP parameter g] 


S6 low-mass 


VSR2-3 low-mass 


representative spectra 


sample rate (1/At) 


4096 Hz 


4096 Hz 


16384 Hz 


data block duration (Tbiock) 


2048s 


2048s 


2048s 


number of data segments (Ns) 


15 


15 


1023 


data segment duration (T) 


256s 


256s 


4s 


stride (A) 


524288 


524288 


32768 



p] gives a list of parameter names and symbols, which are the same names and symbols used in g] . The second 
and third columns gives the values of these parameters used in S6/VSR2-3 CBC searches. These parameters can be 
used to reproduce the inspiral horizon distance data accompanying this note. The fourth column gives the values 
of the parameters used to compute the representative spectral density curves shown here. In making our choice of 
parameters for computing representative spectra, we sacrificed frequency resolution (A/ = 1/T) for PSD accuracy 
(which increases with Ns). 

One potential pitfall with using these spectra is that the choice of representative PSD for a detector is not obvious. 
Here we illustrate the degree to which our choice of using the mode affected the chosen PSD. We compare the spectra 
for HI corresponding to times when HI operated near its mode to times when it operated near its mean and maximum 
of its inspiral horizon distance distribution. In Table \H\ we provide a quantitative summary of the low-mass inspiral 
horizon distance distributions for a 1.4-1.4M0 binary given in Mpc. We see that the horizon distance varies by roughly 
10% between its mode and mean. This suggests that spectral density curves from a detector's most common sensitivity 
(mode) may differ significantly from the spectral density of a detector's "average" performance. To illustrate this 
point, we plot in Fig. [4] three spectra for HI from different times in S6. 

All of the data used here have been computed using the final version of calibration used in the CBC searches. Note 
that the noise spectra presented here are subject to systematic uncertainties associated with the strain calibration. 
These uncertainties can be up to ±15% in amplitude. For more detail, see references [TOllllj. 



Representative Spectra for LIGO/Virgo Detectors in S6/VSR2-3 
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FIG. 3: Representative spectral density curves for LIGO and Virgo detectors during S6 and VSR2-3. We plot here the 
amplitude spectral density, which is the square root of the power spectral density, since the strength of a gravitational wave 
signal is proportional to the strain induced in the interferometer and the sensitivty of a detector is therefore most naturally 
expressed in terms of this quantity. These spectral density curves correspond to May 9, 2010 (GPS 957474700) for HI, February 
27, 2010 (GPS 951280082) for LI and August 30, 2009 (GPS 935662133) for VI. These times are chosen such that the inspiral 
horizon distance for each detector at that time coincides with the mode of its inspiral horizon distance distribution, as given 
by the midpoint of the most populated bin in Fig. [TV). 

TABLE II: Inspiral Horizon Distance Summary for a 1.4— 1.4-Mq Binary 





HI LI 


VI 


mean 


39.5 34.0 


16.6 


max 


49.3 47.2 


23.2 


mode 


44.1 36.6 


18.8 


std 


4.7 6.9 


3.9 



Acknowledgements 



The authors gratefully acknowledge the support of the United States National Science Foundation for the construc- 
tion and operation of the LIGO Laboratory, the Science and Technology Facilities Council of the United Kingdom, the 
Max-Planck-Society, and the State of Niedersachsen/ Germany for support of the construction and operation of the 
GEO600 detector, and the Italian Istituto Nazionale di Fisica Nucleare and the French Centre National de la Recherche 
Scientique for the construction and operation of the Virgo detector. The authors also gratefully acknowledge the sup- 
port of the research by these agencies and by the Australian Research Council, the Council of Scientic and Industrial 
Research of India, the Istituto Nazionale di Fisica Nucleare of Italy, the Spanish Ministerio de Educacion y Ciencia, 
the Conselleria dEconomia Hisenda i Innovacio of the Govern de les Illes Balears, the Foundation for Fundamental 
Research on Matter supported by the Netherlands Organisation for Scientic Research, the Royal Society, the Scottish 
Funding Council, the Polish Ministry of Science and Higher Education, the FOCUS Programme of Foundation for 
Polish Science, the Scottish Universities Physics Alliance, The National Aeronautics and Space Administration, the 
Carnegie Trust, the Leverhulme Trust, the David and Lucile Packard Foundation, the Research Corporation, and the 



10 



10 



Representative Spectra For HI During S6 



tq 



10 ■■ 



ra 10 ■' 



i 

Q- 
IJ1 
(LI 



|- 10" 22 
< 



io- 





■ ■ 1 














, , , 1 






, , 


: 






— mode ■ 

— mean ; 
best 




























I 






w. 








1 


w 


yi 


P 


. 


I 


u 


£ 


J^* : 


Nisi] 










>* 





10 2 



10 J 



Frequency (Hz) 



FIG. 4: Various possibilities for a representative PSD for HI during S6. These spectral density curves correspond to times 
when the detector operated near its S6 mode (44.1Mpc), mean (39.5Mpc) and maximum (49.3Mpc) inspiral horizon distance. 
The times chosen are May 9, 2010 (GPS 957474700) for the mode, November 4, 2010 (GPS 941365351) for the mean and July 
4, 2010 (GPS 962268343) for the maximum. 



Alfred P. Sloan Foundation. This document has been assigned LIGO Document No. ligo-t 1100338. 

Anyone using the information in this document and associated material (S6/VSR2/VSR3 noise spectra, inspiral 
ranges, observation times) in a publication or talk must acknowledge the US National Science Foundation, the LIGO 
Scientific Collaboration, and the Virgo Collaboration. Data files associated with the results and plots presented in 
this document can be found here: https://dcc.ligo. org/cgi-bin/public/DocDB/ShowDocumcnt?docid=63432. Please 
direct all questions to the corresponding author (Stephen Privitera, stephen.privitera@ligo.org). Please inform the 
corresponding author and the LSC and Virgo spokespeople (LSC-Spokesperson@ligo.org and Jean- Yves. Vinet@oca.eu, 
respectively) if you intend to use this information in a publication. 



11 



[1] The LIGO Scientific Collaboration and Virgo Collaboration: B. Abbott, ct. al. Sensitivity Achieved by the 

LIGO and Virgo Gravitational Wave Detectors during LIGO's Sixth and Virgo's Second and Third Science Runs. 

http://arxiv.org/abs/1003.2481. 
[2] LIGO Scientific Collaboration and Virgo Collaboration: J. Abadie, J., et. al. Search for Gravitational Waves from Low 

Mass Compact Binary Coalescence in LIGO's Sixth Science Run and Virgo's Science Runs 2 and 3. Phys. Rev. D., 2012. 
[3] See lines 964-990 in the lalapps tmpltbank source code at 

http://www.lsc-group.phys.uwm.edU/cgit/lalsuite/tree/lalapps/src/inspiral/tmpltbank.c 

and the function compute_candle_distance in lines 86-119 from 

http://www.lsc-group.phys.uwm.edU/cgit/lalsuite/tree/lalapps/src/inspiral/inspiralutils.c. 
[4] Bruce Allen, Warren G. Anderson, Patrick R. Brady, Duncan A. Brown, and Jolien D. E. Crcighton. FINDCHIRP: an 

algorithm for detection of gravitational waves from inspiraling compact binaries. arXiv:gr-qc/0509116. 
[5] Lee Samuel Finn and David F. Chernoff. Observing binary inspiral in gravitational radiation: One interferometer. Phys. 

Rev. D, 47:2198-2219, 1993. 
[6] Duncan A. Brown. Searching for Gravitational Radiation from Binary Black Hole MACHOs in the Galactic Halo. 

http://arxiv.org/abs/0705.1514. 
[7] Patrick J. Sutton. S3 performance of the ligo interferometers as measured by sensemonitor. LIGO-T030276-vl, 

https://dcc.ligo.org/cgi-bin/private/DocDB/ShowDocument?docid=27267. 
[8] LIGO and Virgo Collaborations. Sensitivity to Gravitational Waves from Compact Binary Coalescences 

Achieved during LIGO's Fifth and Virgo's First Science Run. T0900499-vl9, https://dcc.ligo.org/cgi- 

bin/DocDB/ShowDocument?docid=6314. 
[9] The LIGO Scientific Collaboration and Virgo Collaboration: J. Abadie, et. al. Search for gravitational waves from binary 

black hole inspiral, merger and ringdown. Phys. Rev. D., 83:122005, 2011. 
[10] S6 Calibration Homepage. https://wiki.ligo.org/DAC/S6Calibration. 
[11] The Virgo Collaboration: T. Accadia, et. al. Class. Quantum Grav., 28:025005, 2011. 



